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Abbreviated title 



ASYMPTOTIC APPROXIMATIONS FOR ELLIPTIC INTEGRALS 



1 Introduction 



A recent table of elliptic integrals p|, |10|, [TT|, [12], uses symmetric standard integrals 
instead of Legendre's integrals because permutation symmetry makes it possible to unify 
many of the formulas in previous tables. Fortran codes for numerical computation of 
the symmetric integrals, which are homogeneous functions of three or four variables, can 
be found in several major software libraries as well as in the supplements to @, |10| 



For analytical purposes it is desirable to know how the homogeneous functions behave 
when some of the variables are much larger than the others. For all such cases we list in 
Section 2 asymptotic approximations (sometimes two or three approximations of different 
accuracy), always with error bounds. Proofs are discussed in Section 3. In most cases 
the approximations are obtained by replacing the integrand by a uniform approximation. 
Many of the results found by a different method in ||16[ have been improved by sharpening 
the error bounds or by finding bounds for incomplete elliptic integrals that are still useful 
for the complete integrals, which are then not listed separately. Cases not considered 
1^] include two for a completely symmetric integral of the second kind and two for 



m 



a symmetric integral of the third kind in which two variables are much larger than the 
other two. 

We assume that x, y, z are nonnegative and at most one of them is 0. The symmetric 
integral of the first kind, 

1 

RF{x,y,z) = - [{t + x){t + y){t + z)]-'/^dt, (1) 
2 Jo 

is homogeneous of degree —1/2 in x, y, z and satisfies Rf{x, x, x) = x"^/^. The symmetric 
integral of the third kind, 

3 

Rj{x,y,z,p) = - [{t + x)it + y){t + z)]-^/\t + p)-'dt, p > 0, (2) 
2 JO 

is homogeneous of degree —3/2 in x,y,z,p and satisfies Rj{x^ x, x, x) = x~^^'^. If p = 
z, Rj reduces to an integral of the second kind, 

RD{x,y,z)=Rj{x,y,z,z) = - [{t + x){t + y)]-^/\t + z)-^/^dt, z > 0, (3) 



2 Jo 



2 



which is symmetric in x and y only. If two variables of Rp are equal, the integral becomes 
an elementary function, 

1 

Rc{x,y) = RF{x,y,y) = - {t + x)-^'\t + yY^dt, y > 0. (4) 

2 JO 

If X < y it is an inverse trigonometric function, 

Rc{x,y) = - arccos(x/?/)^/^, (5) 

and if X > 2/ it is an inverse hyperbolic function. 



R^{x, y) = {x- y)-^/\Tccosh{x/yy/^ = {x - y)-^'^ In ^ - (6) 

If the second argument of Rc is negative, the Cauchy principal value is (4.8)] 

/ \ 1/2 

Rc{x,-y)=(^^] Rc{x + y,y), y > 0. (7) 
\x + yj 

If the fourth argument of Rj is negative, the Cauchy principal value is given by |T8| , (4.6)] 

{y + p)Rj{x,y,z,-p) = {q - y)Rj{x,y, z, q) - 3RFix,y, z) 

/ xyz \ 

+ 3 ■ Rc{xz + pq,pq), p > 0, (8) 

\xz + pq J 

where q — y = {z — y){y — x)/{y + p). If we permute the values of x,y,z so that 
{z — y){y — x) > 0, then q > y > 0. 

A completely symmetric integral of the second kind is not as convenient a.s Rd for use 
in tables because its representation by a single integral is more complicated |^, (9.1-9)]: 

It is symmetric and homogeneous of degree 1/2 in x,y,z, and it satisfies Rg{x,x,x) = 
x^/^. It has a nice representation by a double integral that expresses the surface area of 
an ellipsoid 0, (9.4-6)]. It is related to Rd and Rp by (|^) and by 

6RG{x,y,z) = x{y + z)RD{y,z,x) + y{z + x)RD{z,x,y) + z{x + y)RD{x,y,z), (10) 
6RGix,y,0) = xy[RDiO,x,y) + RniO,y,x)]. (11) 



Legendre's complete elliptic integrals K and E are given by 

K{k) = RpiOA-k^l), (12) 
E{k) = 2i?G(0, 1 - A;^ 1) 

= ^^^[RDiO,l-k^l) + RD{OAA-k% (13) 

K{k)-E{k) = yi?^(0,l-fc2,l), (14) 
E{k)-{l-e)K{k) = :^!il-l^i?^(0,l,l-fc2). (15) 

Approximations and inequalities for K, E, and some combinations thereof are given in 
|l|, H, 1^. If the error terms in (pOf), (^TJ), and (^) are omitted, the approximations reduce 
to the leading terms of well-known series expansions of K and E for k near 1 p!5| , p. 54] 
0, 900.06, 900.10]. If the series for K is truncated after any number of terms, simple 
bounds for the relative error are given in [|14|, (1.17)]. A generalization of this series to 
Rf{x, y, z) with x, y « z is given in |1^, (1.14)-(1.16)], again with simple bounds for the 
relative error of truncation. 

The various functions designated by R with a letter subscript are special cases of the 
multivariate hypergeometric i?- function, 

^-a(fel,..., "n ) ^1 ; • • • ; ) ; 

which is symmetric in the indices 1, . . . ,n and homogeneous of degree —a in the variables 
zi, . . . , Zn- Best regarded as the Dirichlet average of x'" [0, § 5.9], it is a symmetric 
variant of the function known as Lauricella's F^. By the method of Mellin transforms, 
series expansions are obtained in H, (4.16)-(4.19)] that converge rapidly if some of the 
z^s are much larger than the others and if the parameters satisfy Z^ILi bi > a > 0. Thus 
the leading terms of these series provide asymptotic approximations for all except Rq 
among the functions 

Rpix.y.z) = R_i{l,^,l; x,y,z), Rc{x,y) = R_^{\,1] x,y), (16) 
Rj{x,y,z,p) = i?_|(i,i,i,l; x,?/,z,p), RD{x,y,z) = R_i{l,^, I; x,y,z), (17) 
Rcix.y.z) = /?i(|,i,i; X,?/, z). (18) 



However, error bounds for the approximations are more easily derived by the methods of 
the present paper. Another function that is used repeatedly in obtaining error bounds is 
0, Ex. 9.8-5] 

R.i{\,\,l-x,y,z) = [{t + x){t + y)r^'\t + zr^dt (19) 
= 2Rc{{^+z)\{^+^fz). (20) 

In Section 4 the asymptotic approximations are applied to show that RF{x,y, z), 
RDix,y, z), Rj{x,y, z,p), and (xyz)'^/"^ are linearly independent with respect to coeffi- 
cients that are rational functions of x, y, z, and p. An Appendix contains some elementary 
inequalities that are used in obtaining error bounds. 

The results in this paper provide upper and lower approximations that approach the 
elliptic integrals as selected ratios of the variables approach zero. Approximations that 
approach the integrals as all variables approach a common value have been found by 
other methods. For example, the theory of hypergeometric mean values yields upper and 
lower algebraic approximations for all the integrals in this paper Thm. 2], while trun- 
cation of Taylor series about the arithmetic mean of the variables gives approximations 
with errors that may be positive or negative. Successive applications of the duplication 
theorem for Rp, making its three variables approach equality, provide ascending and de- 
scending sequences of successively sharper (and successively more complicated) algebraic 
approximations to Rp and Rc 0]. Transcendental approximations that approach Rp 
when only two of its variables approach equality are furnished by 

Rc{x,^^^ <RF{x,y,z)<Rc{x,^), yz^O, (21) 

which follows from (|7lD. The inequalities can be sharpened by first using Landen or Gauss 
transformations of Rp [0, § 9.5] to make y and z approach equality. If x = the Gauss 
transformation reduces to replacing ^/y and ^/z by their arithmetic and geometric means, 
and each i^c-function becomes 7t/2 divided by the square root of its second argument. 
Therefore, in the complete case the procedure reduces to the algorithm of the arithmetic- 
geometric mean 0, (6. 10-6) (9.2-3)] and provides ascending and descending sequences of 
algebraic approximations, of which leading members are shown in (|33|). 

5 



2 Results 



We assume throughout that x, ?/, and z are nonnegative and at most one of them is 0. 
The last argument of Rq-, Rd-, and Rj is assumed to be positive (see and (||)). 



CI. Rc{x,y) with x « y. 



Rcix,y) 



TT 



TTX6 



2^ y %3/2 



where 1/(1 + Jx/y) < 9 <1 with equahties iff a; = 0. 



(22) 



C2. Rc{x, y) with y « a;. Two approximations of different accuracy are 

'=fln^ + :^ln^) 



Rc{x,y) 



2v/i 
1 



y 2x-y 



y 



Sy' 



, , , 4x y 
1 + — In —+ , 

2xJ y 2x Ax[2x — y) 



In 



Ojx 

y 



where 1 < 0j < 4 for i = 1,2. The first approximation imphes 



Rcix,y) < 



1 



In 



Ax 



2y^{l-y/2x) y- 



(23) 
(24) 

(25) 



Fl. Rf{x, y, z) with £C, y « z. Let a = {x+y)/2 , (7 = y^x^ , and p = max{x, y}/z. 
Then 

i?ir(x,?/, Z) 



1 



r \ 

In + 7^ h 



where 



^ 1 2^ ^ 
m < r < 



l-g/z a + g 
The upper bound imphes 

RF{x,y,z) < 



a , oz 
in 



1 — a/2z a + (? 



In 



8z 



2^{\- al2z) a + g' 
A sharper lower bound and a higher-order approximation are given by 

1 / . 82 ari 



RF[x,y,z) 



r- (in + 

2^ \ a + g 2z 

1 



2^z 



^ ^ a \ 82; 2a — g ^ 3(3a^ — 5'^)r2 



22/ a + 



2z 



16z^ 



(26) 



(27) 

(28) 
(29) 
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where 



, ^ ln(l/p) 1 , Sz 

In — < ^Irl <ri< — In , 2 = 1,2. 

2a 1 — p l — a/2z a + g 



By (|T2|) this imphes (since 4^^ < 4 - /c'^ if A;^ < 1) 



/ 4 9k'^ O2 



where < fc' = Vl - and 1 < < 4 for z = 1, 2. 



F2. RF{x,y, z) with z « x,y. Let a = {x + y) / 2 and g = ^/xy . Then 



Rf{x, y, z) = Rf{x, y,0)-^ + (32) 

where l/{l + ^fz/g) <e < a/g. Note that Rp^x.y.Q) = n/2AGM{y/x, ^y), where AGM 
denotes Gauss's arithmetic-geometric mean 0, (6. 10-6) (9.2-3)], and hence 

1 I 2 2 2 / 2 \ 1 

Va + 9 ^[a + g)/2 + ^ \ag + g^ y/g 

with equahties iff a; = y. 



Dl. Ru{x , y , z) with x,y <.< z. Let a = (x -|- ?/)/2 and (7 = y^a^. Then 



where 



Rn{x, y,z) = -^(ln^-2 + -ln , (34) 

2z-^'^ \ a + g z a + g 



<e < 



l-g/z 2{l-a/z) 



D2. R£){x , y , z) with z<.<Cx,y. Let a = (x + ?/)/2 and (7 = ^/a^. Then 

,/xyz \ 2 y ^ 



RD{x,y,z) = ^(l-^,^], (35) 
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where 

1 4 a 

1- -J- <e < -. 

V ^ 9 

A higher-order approximation is 

Rd{x, y, z) = Rd{0, X, y) - Rd{0, y, x) + 



2g\l + ^z/g) 
where 

1 „ 3a 
<d < 



/2/3 + ^z/g 2g{l + ^z/g) 

An approximation of still higher order is 



Roix, y, z) = —= Rg{x, y, + \\- - 

^xyz xy g-^ \ A \ a 



where we have used (jO) and where 



1 + Jz/a \9, 



9 I ' 



D3. R£){x,y, z) with y,z<.<Cx. Let a = + 2;)/2 and (7 = y^y^. Then 

RD{x,y,z) = [— ^ ) , 

^yx \g + z Ax J 

where 

1 ^ 2x 2z 1 ^ 8x 
In < r < ; — In 



l — g/x a + g g + z l — a/2x a + g 



D4. Ruix,y, z) with x«y,z. Let a = + 2;)/2 and (7 = Then 

Rd{x, y, z) = RniO, y,z) + -^[-l + —J- 

gz \ 4 V a 



5 



where 



1 + Jx/a \9, 



a 



Jl. Rj(x,y, z,p) with x,y,z«p. Let a = {x + y + z) /3 and b = {^/3/2){xy + 
xz + yzY^"^ . Then 

3 Svr 
Rj{x,y,z,p) = -RF{x,y,z) + ^-^(-1 +r), (40) 

where 



h/p 3 ^a/p 

< r < 



1 + ^b/p 2 1 + ^a/p 

In the complete case a sharper resuh is 

RAx, y, 0. p) = 5 (^R,(,, y, 0) - (i + ^) . (41) 

where y/xy < < {x + y)/2 with equahties iff x = 



J2. Rj{x,y, z,p) with p « x,y,z. Let (? = (xt/z)-*^/^, 3/;."-'^ = x^^ + y^^ + z~^, 
and A = ^^xy + ^/xz + Note that g is the geometric mean and h is the harmonic 

mean, whence g > h with equahty iS x = y = z. Then 

i?j(x, z, p) = ( In ^ - 2 + r ) , (42) 



where 

h 2{g -p) p 

A higher-order approximation is 



m — < r < — m — . 



where 



i?j(x, y, z,p) = In ^ + 2i?j(x + A, y + A, 2; + A, A) + (43) 

2y/xyz py A^xyz 



2 n 3 , /i 
m — < r < m — . 



g — p p h — p p 

The second term in the approximation is independent of p but is otherwise as compli- 
cated as the function being approximated. The same is true of an even more accurate 



approximation |16, Thm. 11] in which the error is of order p instead of p\np and the 
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leading term involves Rc- 



33. Rj(x, y, z, p) with x, y « z, p. Let a = {x + y)/2 and g = Jxy. Then 



2^p 



In- 



a + g 



p\ 6 2p 



-2Rc h-] + -ln 



z ) p a + g 



(44) 



where 



9 n 

<e < 



^~9/p l — a/p\ 2zJ 



1 + 



p 



34. Rj{x, y, z,p) with z, p « x, y. Let a = {x+y)/2, g = ^/xy, b = \j3p{p + 2z)/2, 
and d = {z + 2p)/3. Then 



Rj{x,y,z,p) = -Rc{z,p) ^ 

9 9-P 



p 

Rc{z,g) - -Rc{z,p) 
9 



(45) 



where I < < a/ g with equalities iS x = y. Since z « g, Rc{z,g) can be estimated 
from (|22|). In the complete case (H3) reduces to 



Rj{x,y,0,p) 



3n 



1 



2^fxyp\ v^+VP/ 
with 9 as before. A higher-order approximation is 

Rj{x, y, z, p) = -Rciz, p) Rg{.x, y, 0) + - — , 

9 2x?/ 



(46) 



(47) 



where we have used (|Tl]) and where 

Vb 



<6 < 



3a ^/d 



1 + Jb/g 2(7 1 + h/g 



35. Rj(x, y, z, p) with x « y, z, p. Let a = {y + z)/2 and g = Jyz. Then 



Rj{x, y, z, p) = Rj{0, y, z, p) + 



V ^ 9) 



(48) 



10 



where 

1 + ^x/a 9 P 



36. Rj{x,y, z,p) with y,z,p«x. Let a = {y + z)/ 2 and g = ^/yz. Then 



Rj{x,y,z,p) = 

'X 



Rciig + pY,2{a + g)p)-- 



where 



1 , 2x 2p ^ , , .0 , N N 1 , 8x 
In -Rc{{9 + p) , 2(a + g)p) <r< hi ■ 



X — g a + g x 
In the complete case this reduces to 



X — a/2 a + g 



3 3s 
Rj{x, 0, z, p) = ^=Rcip, z) - 



xp 



4x3/2' 



where 



In 



Ax 



— -2^Rc{p,z) <s<Y 



In 



16a; 



z/Ax z 



(49) 



(50) 



Gl. Rg{x, y, z) with x, y « z. Let a = {x + y)/2 and g = Jry. Then 



i?c(x,2/,.) = ^(l + ^), (51) 
where 

a + g , 2z 4a , , , 2z a 

— ^In^— + 2(7- — <r < 3a-(7 \n^—^2g--. 
2 a + g 3 a + g 3 

In the right-hand inequality it is assumed that 5a < A sharper result for the complete 
case is 



where 

3, z 1/162 13\ 
- In - < s < — In . 

Ay l-y/zy y 6 J 

By (|13D this follows from 

i?(fc) = l + ^(lni-l + fc'V), (53) 
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where < A;' = \/l — k'^ « 1 and 



3, 1 1 /, 4 13\ 

— m — < r < In 

8 k' k{l + k)\ k' 12 J 



G2. Rg{x, y, z) with z « x, y. Let a — {x + y)/2 and g — Then 

Rg{x, y, z) = Rg{x, y, 0) + 7i9z/8, (54) 

where 

1 /A ^\ / 2 V^' 1 



11--J- <9<[^-^ < 
Va V TrVa/ \ag + gy y/g 

with equahty iff x = y. 

3 Proofs 

Most of the results in Section 2 are obtained by replacing an integrand / by an ap- 
proximation fa, writing / / = I fa + f{f ~ fa), and finding upper and lower bounds for 
Kf ~ fa)- AH integrals are taken over the positive real line. The function fa is usually 
chosen to be a uniform approximation fa = fi + fo — fm, where fi is an approximation in 
the inner region, in the outer region, and in the overlap region or matching region. 
For instance, if f{t) — [{t + x){t + y){t + z)]~^/'^ with x, y « z, we get fi by neglecting 
t compared to z, fo by neglecting x and y compared to t, and by doing both. A first 
example of this process is the proof of Lemma 1. 

Lemma 1. If x > 0, y > 0, and < x + y « z, then 



f 

Jo 



* ' In^, (55) 



'{t + x){t + y){t + z) « + ^' 

where ^/xy — g < 9 < a — {x + y) /2 with equalities iff x — y 
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Proof. Let 

fit) = , ^ — , m= , ^ 

^it + x){t + y){t + z) z^{t + x){t + y) 

1 1 

foif) = r, fmit) = ~ ■ 

t{t + Z) zt 

Taking fa = fi + fo-fm, we find 

1 , 2z 



/ fa{t)dt = -In 

Jo z 



a + g 



and 

til 1 



/ /a 



Inequality (|6^ in the Appendix implies 

Q 

z^{t + x){t + y){t + z) 

and thus 

f=ffa+[{f-fa)=f fa + - I f=-r\r I ^'^ = ^1^^- ° 

J J J zJ 1 — d/zj z — d a + g 

As a second example, in which Lemma 1 is used, consider Rf{x, y, z) with x, y « z. 

Let 

fit) = , m ^ 



^{t + x){t + y){t + z) ^{t + x){t + y)z 



at) = —l=, ut) ^ 



t^/t + z' ™ 
Taking fa = fi + fo — fm, we find (with a and the same as before) 



roo I 

/ /,(t)rft = ^ln 
JO \/z 



\fz a + g 

and 

/ fa 



Inequalities (pi]) and imply 

</-/.< 



2^^(t + x)(t + 2/)(t + z) 22^(t + a;)(t + 2/)(t + z)' 
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Hence, by Lemma 1, 
9 



In 



2z 



fa 



2y^{z — g) a + g J 2z J 1 — a/2z 

where the last inequahty follows from the next to last. We complete the proof of (^) by 
noting that 

2RF{x,y,z) = I f = I + l(f-f^). 

Equations (^) and (|29| ) are obtained from |T^, (2. 15) (3.25)] with w = oo. To derive 
(^) we construct fa = fi + fo — fm as usual and find bounds for /(/ — fa) by using (|60|) 



and (|65D . To simplify the upper bound we note that RF{x,y,z) < RF{x,y,0) and use 



Equations (H) , (H) , (|2|) , and (||) follow from {^,(^,(^, and respectively, 
by replacing x by y, replacing z by x, and simplifying. 

Among the approximations for Rd we need discuss only ( pSj) and (p7|), since (|3^), 
and (|39|) follow from (|4^) , (l47|), (|49|) , and (|48|) , respectively, by putting p = ^ 



' 5 VMM/ 5 



and simplifying. To prove (^) we let 
fit) ^ 



^{t + x){t + y){t + z)^/^ 
choose fa = fi, and apply (|^) to get 

t 



f^{t) 



1 



1 



g{t + g){t + z)^^ 
(9 z \ 



g{g - z)^/t + z \t + g t + 



Rciz,g) 



< fa-f< 



< fa-f< 



g{t + z)3/2 



at 



g2j{t + x){t + y){t + z)^l^ 



9- z 



9 



g^Vt + z{t + g) 
2a 
? 



r 2a 

< J{fa-f)<-Rc{z,9). 



Use of (p2D completes the proof. Approximation (37) follows from applying (^) to two 
terms on the right side of 



Rd{x, y, z) = 3{xyz) - Rd{z, x, y) - Rd{z, y, x) 



(56) 



an identity that comes from 0, (5.9-5) (6.8-15)]. 



14 



In discussing approximations for we define 

fit) = , ^ 

^{t + x){t + y){t + z){t + p) 

and construct fi, fo, and fm for each case in the manner described at the beginning of 

this Section. For example, if x, y, z « p, then fi is obtained by neglecting t compared 

to p. Unless otherwise stated, we define fa = fi + fo — fm, take / fa as an approximation 

to / /, and find bounds for J{f — fa) by using the inequalities in the Appendix. 

To prove pO|) we use ([691). To prove (0) we use (Q) and note that /(/ — /a) = 



{6/p) J f. Before discussing (^2]), we consider (|4^), in which the error bounds are easily 
found by using (|7^). Finding / fa requires an integration by parts and a formula of which 
we omit the proof, 

/ {lnt)-[{t + x){t + y){t + z)]-'/^dt = ^=\n- Rj{x + X,y + X, z + X, X), (57) 

JO at y/xyz Axyz 6 

where A = ^Jxy + ^/xz + ^Jyz. To have a simpler approximation (^), we define fa = 

fi + fo - fm and 4>a = fi + fs - fm, wherc fo has been replaced by 

Then 



y/xyz \ p 



and an upper bound for J{f — (pa) is found by using ^ (t + x){t + y) (t + z) > {t + g)^/"^ 
and (0). To find a lower bound, we note that f — fa > 0, whence 

f - (Pa = f - fa + fo- fs> fo- fs- 

A lower bound for /(/o — fs) follows from (ffSl) . 

The straightforward proof of (^4]) uses (|64|), (|67D, and Lemma 1. For the elementary 
approximation ( ^5D we choose fa = fi and use For the more accurate approximation 
(^) we take fa = fi + fo — fm and evaluate / fa by integrating by parts. The error bounds 
follow from ( |66D and ( p9| ) with two variables equated. To find the error bounds for P8|), 
we use (pDD, and ( [7TD to prove 



(t + a;)(t + a) x \g pj y/t + z{t + g) 

15 



After integrating, (^) is used to complete the proof. In the case of (|^), where J{fo — fm) 
is infinite, we choose /„ = fi and evaluate / fa by (^). It follows from (^) that 

^ <fa-f < 



2^x{t + y){t + z) V + ^ x{t+p)J 2x^{t + x){t + y){t + z)' 

where we have replaced t/{t + p) by 1 in the upper bound and x/(a; — p) by 1 in the lower 
bound. We then use (||), and (^. 

The function Rq can be expressed in terms of Rp and Rd by (|17D and Table 9.3-1]: 

1 j xy 

2RG{x,y,z) = zRF{x,y,z) - -{z - x){z - y)RD{x,y, z) + J — . (58) 

Applying (^) and (p^), we obtain (pT]). The error bounds have been substantially sim- 
plified by using the numerical value of In 2 and assuming 5a < 2; in the upper bound. It 



is not hard to obtain ( ^3]) from a well-known infinite series |15, p. 54] for E{k) by using 
the inequality 

1 + h'"" < 2Fia, |; 2; k'') < (1 - k'')-^/' = 1/k, < k' < 1, 
8 

for the hypergeometric function 2-^1- Unfortunately (^) does not lead to simple error 
bounds for (|5^). Instead, we define f{z) = RG{x,y,z) and find from 0, (5.9-9)(6.8-6)] 
that 

f [z) = 



'0 ^{t + x){t + y){t + zfl^' 
Since this is a strictly decreasing function of 2;, the mean value theorem yields f{z) 
/(O) + zf\0 where 



/'(^)<r(C)<r(0) = ^i?Hx,2/,0). 



By (|7TD and (|) we see that 



1 /""^ tdt 1 



Use of (^) and ( p2D completes the proof of (p^. 
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4 Application to linear independence 

In 0, Thm. 9.2-1] it is shown that RF{x,y, z), RG{x,y, z), an integral of the third kind 
called RH{x,y, z,p), and the algebraic function (xyz)^^^'^ are linearly independent with 
respect to coefficients that are rational functions of x,y,z,p. It then follows [j^, §9.2] 
that every elliptic integral can be expressed in terms of Rp, Rq, Rh, and elementary 
functions. From (^) and a known relation expressing Rh in terms of Rj and Rp, we 
may conclude that every elliptic integral can be expressed in terms of Rp, Rd, Rj, and 
elementary functions. In order to reach the same conclusion without invoking Rq and 
Rh, we shall use the results of this paper to prove the linear independence of Rp, Rd, Rj, 
and {xyz)~^^'^ with respect to coefficients that are rational functions. 

Theorem 1. The functions Rp{x,y, z) , R£){x,y,z), Rj{x,y,z,p), and {xyz)~^^'^ are 
linearly independent with respect to coefficients that are rational functions ofx,y,z, and 
P- 

Proof. Let a,/?, 7, and 6 be rational functions of x,y,z, and p. We need to prove 
that 

aRp{x, y, z) + PRd{x, y, z) + 7i?j(x, z, p) + 5 {xyzY^^'^ = (59) 

iff a, 7, and 6 are identically 0. We may assume that these coefficients are polynomials 
since we can multiply all terms by the denominator of any rational function. As p — > 0, 
(^) shows that Rj{x, y, z,p) involves Inp while all other quantities are polynomials in p, 
whence 7 = 0. As z ^ 00 we have 

a = a;z™(l + 0(1/ z)), (3 = ^^"(l + 0{l/z)), 

where m and n are nonnegative integers and a and b are polynomials in x, y, and p. Using 
(^) and (|3^) and multiplying all terms by 2z^^'^, we find 

+ 25{xy)-^/^z = 0. 



az 



m+l 



In 



a + g 



+ 



Inz 



+ 3bz' 



In- 



a + g 



2 + 



Inz 
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Cancellation of the leading terms in Inz requires az"^~^^ + Sbz"^ = 0, implying n = m + 1 
and a = —3b and leaving 

0(z™ In z) - 6bz"'+^ + 26 {xyY^'^z = 0. 

Because the second term is of different order from the first and does not have a square 
root in common with the third, it follows that 6 = 0, whence also a = 0. Since the 
leading terms of the polynomials a and (3 are identically 0, so too are a and (3. Finally, 
with only one term remaining in (p9|), we have 5 = 0. □ 

It is an open question whether Theorem 1 is still true if the coefficients are algebraic 
functions instead of rational functions. However, polynomial coefficients suffice (see the 
first paragraph of § 9.2]) to prove that every elliptic integral can be expressed in terms 
of Rf, Rd, Rj, and elementary functions. 



Appendix 



Elementary inequalities 

Assuming x, y, z, and t are positive, we list and prove some inequalities that are used in 
this paper to obtain error bounds: 

X 1 1 X , , 

< ^ - < (60) 



2Vi{t + x) Vi VtTx 2Wt + x' 
tilt 

< ^ - -^r= < o (61) 



2^/x{t + x) ^Jx y/t + X 2x\/t + X ' 
1 1 _ Ox ^ R ^ 

1<«<|. (63) 



a;3/2 (t + x)3/2 a;3/2(t + x)' 2' 
In the next five inequalities let a = {x + y)/2 and g = ^Jxy. Inequalities become equalities 
in (|e|), (H), and (|§) iflF x = 

, = , =, 9<e<a, (64) 

t ^{t + x){t + y) t^{t + x){t + y) 



til at , , 

< ^ - ; < ; (65) 



g{t + g) ^ ^{t + x){t + y) ^ {t + x){t + y) 
or alternatively, 

' ' i<e<^, (66) 



J{t + x){t + y) 9it + g) 9 



' ' l<^<l + f, (67) 



y/xy ^t + x(t + y) ^Jxyit^y)' 2x' 

1 Qt ^ 0. Q , ^ 

l<e <- + -. (68) 



^{t + x){t + y){t + z) gz^{t + x){t + y)' 9 ^ 

Finally we have 

6 1 1 3a , , 

< TiiT^ - ^— — — < T^:^T7^7^-:T' (69) 



^^/^(^ + &) t^^' ^{t + x){t + y){t + z) 2ty^{t + ay 



where a = {x + y + z)/3 and b = y3{xy + xz + yz)/2, and 

t 1 1 3t 

< ^= - I < n ^l.u , (70) 



^3/2(t + ^) y'(t + x)(t + y)(t + 2) 2g^/^{t + h) 

where g = (xyzY^^ and 3h^^ = x^^ + y^^ + z^^. 
To prove (|60|) we write 

1 1 VtTa; - y/t X 



Vt VtT^ ^t{t + x) ^t{t + x)iVtT^+Vi) 
and replace the last denominator factor by either or 2y/t + x. Interchange of t and 



X 



leads from (|60|) to (|61|) . To prove (|62|) let y = ^1 + x/t and write 



which increases from 1 to 3/2 as t increases from to oo and y decreases from oo to 1. 



Interchange of t and x leads from (|62D to (|63D. 



If the left side of (|6^) is put over a common denominator, it suffices to observe that 



t + g < Jit + x){t + y) <t + a. (71) 
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The left inequality is enough to prove the left inequality in ([651). To prove the right 
inequality in we define 



and note that (f){t) tends to 0/(7 as t —> and to 1 as t ^ 00. Differentiation shows that 
(j) decreases monotonically, because 



(t + + y)^> = ^^ta + + [{ta + g^f - - g")]"^ < 0, 

with equality iff x = Because of (|7TD, ( |65D implies (|6^). 

Equation (pT]) is proved by solving for 6^ and using (pT]). Likewise, (^) is proved by 
solving for 

^ = m 



t + z 



and using the result just established that 1 < (f){t) < a/g. 

To prove (|69D we use Maclaurin's inequality [0, Thm. 52] to find that 



r + 2br + 4bH/3 <{t + x){t + y){t + z) < {t + a) 



and hence 



(72) 



Vi{t + b) < ^ {t + x){t + y){t + z) <{t + af^. 

Inequality (|69D follows from this and (|6^) . 

The proof of ( |70D uses Maclaurin's inequality and the inequality of arithmetic and 
geometric means to get 



t + g 



< 



{t + x){t + y){t + z) 



xyz 



1/3 



~^ x) ( ~^ ( z 



1/3 



< 1 + 



with equalities iS x = y = z, whence 

{t + gf'^ < ^{t + x){t + y){t + z) < (I) {t + hfl\ 
Two applications of ( p^ ) complete the proof of (fTop. 



(73) 
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